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Abstract 

The results from direct numerical simulations of a spatially evolving, super- 
sonic, flat-plate turbulent boundary-layer flow, with free-stream Mach number 
of 2.25 are presented. The simulated flow held extends from a transition region, 
initiated by wall suction and blowing near the inflow boundary, into the fully 
turbulent regime. Distributions of mean and turbulent flow quantities are ob- 
tained and an analysis of these quantities is performed at a downstream station 
corresponding to Re x = 5.548 X 10 6 based on distance from the leading edge. 

1. Introduction 

Data from direct numerical simulations of turbulent shear flows are becoming an 
essential guide to the development of turbulent correlation closure models and the 
analysis of flow dynamics. In the incompressible regime, such simulations of turbu- 
lent channel (Mansour, Kim and Moin 1988) and boundary-layer (Moin et al. 1990) 
flows have been able to identify the relative distribution of the various terms in the 
Reynolds stress and dissipation rate budgets. In addition, analyses of the turbulent 
structures (Moser 1988. Moin and Moser 1989) in turbulent channel flow have also 
been made. These studies have significantly contributed to the understanding of 
turbulence physics as well as the development of improved turbulence models. 

There has been and continues to be interest in the computation of turbulent 
compressible flows motivated by the optimal design of high-speed vehicles and asso- 
ciated propulsion systems. Databases of turbulent flow's, which are more commonly 
available in the incompressible regime, have not been as wudespread for supersonic 
compressible flows; however, numerical simulations of such supersonic flows have 
been performed (Rai, Gatski and Erlebacher 1995, Coleman, Kim and Moser 1995, 
Guarini et al. 2000, and Maeder, Adams and Kleiser 2001). 

The direct numerical simulation (DNS) presented here is of a turbulent, super- 
sonic, zero pressure gradient, flat plate boundary layer flow with a free-stream Mach 
number = 2.25 and a free-stream Reynolds number Reo ~ 4.25 X 10 3 . Unlike 
previous boundary-layer studies (Coleman, Kim and Moser 1995, Guarini et al. 2000, 
and Maeder, Adams and Kleiser 2001), the numerical simulation is a fully spatially 
evolving boundary-layer flow with an inflow boundary located in the laminar regime, 
and transition is forced through wall suction and blowing. The conditions closely ap- 
proximate the experiment of Shutts, Hartwig and VVeiler (1977), Case 55010501. A 
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flow with this large a value of free-stream Reg can be simulated because of the depen- 
dence of viscosity on temperature. Under the adiabatic conditions of the experiment, 
the temperature increases as the wall is approached and the Reg based on local con- 
ditions is correspondingly reduced. Thus, in regions near the wall where resolution 
requirements are the most stringent, the increased viscosity effectively reduces the 
grid density required to compute the flow. 

The computational method used in this study is a high-order accurate, upwind- 
biased, implicit, finite-difference method that was developed (Rai and Moin 1991, 
1993) and used to compute transition to turbulence on a flat plate with a low Mach 
number ( M ^ = 0.1), and subject to high free-stream disturbance levels. The tech- 
nique involves using upwind-biased differences for the convective terms, central dif- 
ferences for the viscous terms, and an iterative-implicit time integration scheme. The 
computational method (Rai and Moin 1991) uses the nonconservative form of the 
governing equations to obtain the solution. Therefore, this approach is confined to 
flows that do not contain discontinuities (shocks). In the present flat plate simu- 
lation, the free stream Mach number is assumed to be sufficiently low to preclude 
eddy-shockiets. The computed results confirm this assumption. 

Results obtained in the fully turbulent region at a distance of x = 8.8 units 
downstream of the leading edge are analyzed. This location corresponds to a Re x 
of 5.548 X 10 6 and an Reg = 1.0 i/ :< , m 4.25 X 10 2 3 (or a wall momentum thickness 
Reynolds number Reg w = U^0jv w ~ 1.343 X 10 3 ). Profiles of various mean veloc- 
ity and thermal variables, as well as turbulent velocity and thermal variables, are 
discussed. 

2. Problem Formulation and Numerical Solution 

The numerical procedure used in the direct simulations presented here is described 
in detail in Rai, Gatski and Frlebacher (1995). As noted in the previous section, an 
important limitation imposed by using a nonconservative formulation of the governing 
equations is that the method can only be used to compute flows that are free of flow 
discontinuities (such as purely subsonic flows). The nonconservative formulation of 
the governing equations was chosen for this investigation because of the efficiency of 
the approach, the availability of a validated computer program that uses the approach, 
and because the flow' does not exhibit discontinuities. The present method is second- 
order accurate in time and fourth-order accurate in space. 

The computational region of interest is shown in Fig. 1. In the streamwise re- 
direction, it extends from x = 4.0 to x = 15.0 for the turbulence simulation. (Note 
that all spatial distances have been normalized by unit length.) the grid spacing in 
the streamwise direction is sufficiently fine to resolve the flow features of interest up to 
x = 9.0. Thereafter, the grid gradually becomes very coarse so that the disturbances 
in the boundary layer rapidly attenuate and a simple boundary condition can be 
used at the exit boundary. The dimensions of the computational region in the wall 
normal y- and spanwise ^-directions are 3.0 and 0.35, respectively. As Fig. 1 shows, 
the boundaries of the computational region consist of the inlet boundary, the exit 
boundary, the lower surface solid boundary, the upper surface symmetry boundary, 
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Upper surface 
symmetry boundary 



and the two boundaries in the spanwise direction. The boundary conditions used at 
each of these boundaries are as follows. 

Inlet Boundary 

The inlet boundary is comprised mainly of a supersonic flow except for the subsonic 
portion of the boundary layer. The values of all the dependent variables are fixed at 
this boundary, and obtained from a similarity compressible boundary-layer solution). 

Exit Boundary 

The exit boundary is a supersonic exit boundary except for the subsonic portion of 
the boundary layer. Ail the variables are extrapolated from the interior to the exit 
boundary for those exit boundary points where the flow Mach number is greater 
than unity. For boundary points where the streamwise Mach number is subsonic, the 
extrapolation process is implemented for all the variables except the pressure. The 
pressure is set equal to the pressure at the grid point at which the Mach number first 
becomes supersonic. To be accurate, this approximation requires that the distur- 
bances in the boundary layer be attenuated before the exit boundary is reached. As 
described previously, this attenuation is accomplished by using a grid that gradually 
becomes very coarse (geometric stretching) in the streamwise direction as the exit 
boundary is approached. 

Upper Surface and Spanwise Boundaries 

A symmetry boundary condition is used at the upper surface boundary. This bound- 
ary condition corresponds to the one used for a channel flow simulation rather than 
a flat-plate simulation. However, the ratio of the width of the channel to the momen- 
tum thickness at x = 9.0 is about 300. t herefore, the computed solution would be 
very close to a flat plate solution. Additionally, the grid coarsens in the ^-direction 
so that disturbances are attenuated before they reach the upper boundary. 
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A periodic boundary condition is used in the spanwise direction. The width of the 
computational region in the spanwise direction was determined with the help of two- 
point correlations of the velocity components in this direction obtained with coarse 
grids and wider domains. The width L z = 0.35 was chosen such that the auto- 
correlation coefficient values were sufficiently small ((9(10 -2 )) at the tail end of the 
correlation curve. 

Lower Surface Solid Boundary 

The lower boundary of the computational domain corresponds to the plate surface. 
Except for the portion of the plate where blowing and suction are implemented, a 
no-slip boundary condition together and a temperature condition (isothermal wall 
temperature) is specified. The dependent variables are obtained from 

«( i) = v (i) = w (i) = 0, '/(i) = Tun P( l) = P{ 2 ) (1) 

where the subscripts (1) and (2) correspond to the wall and the first point away from 
the wall, respectively. The density at the wall is computed by using the pressure 
and the temperature. In the region of wall blowing and suction (in this computation, 
4.5 < x < 5.0) the normal component of the velocity is computed as 

vi = Au 00 f(x)g{z)h{t) (2) 

= 4 sin 0( 1 — cos 9)1 \/27 
= 27r(a- - x a )/{x b - x a ) 

l max 

= sin( 27 r l(z/z max + (j>i)) 

1=0 

= 1, Z, = 1.25Z l+1 (3) 

TTimax 

= T m srn{2n(pt + 

m = 1 

= 1, T m = 1.257 m+ i 

= 1.5 and xi = 5.5 are the streamwise locations of the 
beginning and the end of the blowing/ suction strip, respectively, A is the amplitude 
of the imposed disturbance (0.04), (3 is the fundamental temporal frequency of the 
disturbance (75000 Hz), <j>i and <b m are random numbers ranging between 0.0 and 
1.0, and z max is the dimension (0.35) of the computational region in the spanwise 
direction. 

file numerical method and boundary conditions just discussed were validated by 
studying the evolution of two-dimensional, small amplitude disturbances in a flat- 
plate flow at a Mach number of 4.5 and a free-stream temperature of 110° R, with 
a unit Reynolds number based on free-stream conditions of 83333. It was shown in 


where 
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0 

M 
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Rai, Gatski and Erlebaclier (1995) that the streamwise evolution of the disturbance 
amplitude maxima computed with the present methodology was in excellent agree- 
ment with parabolized stability equation results. With the success of this small am- 
plitude disturbance study, the methodology was then applied to the forced, flat-plate 
supersonic boundary-layer flow. 

The results reported by Rai, Gatski and Erlebaclier (1995) were obtained in sev- 
eral stages by interpolating an initial solution from an extremely coarse grid onto 
successively finer grids. The computation on each of the intermediate grids was car- 
ried to a point where the skin-friction distribution (after spanwise averaging) showed 
minor changes in time and exhibited a time average that was essentially in equilib- 
rium. In their study, these results were mapped onto the finer grid (N :c X N y X A' ~ ) = 
(1701 X 55 x 501) used in this study and allowed to evolve until the transients from 
the mapping process disappeared. 

To help assess the accuracy of the numerical simulation, one-dimensional spectra of 
the primitive variables were examined in the periodic 2 -direction. At a given x-station 
along the solid boundary and a y-value in the boundary layer, the corresponding two- 
point correlation function can be formed in the homogeneous 2 -direction, and is given 

by 

N;~ 1 

Raai^z) ^ ( &k(Xk-\-k r i i'V 0, 1, . . . , h 1, (4) 

fc=l 

where r z = k r Az , N z is the number of grid points in the 2 -direction, and a represents 
any one of the primitive variables p., u, t\ w. or p. (Note that the periodicity allows 
the two-point correlation function to be obtained over the whole range of 2 grid 
points.) The discrete power spectrum is then obtained from Eq. (4) as 

(AT-D/2 

E aa (k z Az) = 1+ 2 Raa{k r Az) cos 

k r = 1 

where k z Az = n/( N z — 1). and Ax = L~/(N Z — 1). Figure 2 shows the normalized 
spectra for the primitive variables in the homogeneous 2 -direction at three different 
y + locations. (As will be discussed further in Sec. 3., these results are obtained at 
a downstream station corresponding to Re x = 5.548 x 10 6 , or x = 8.8 units from 
the leading edge.) The wall normal distances correspond to an inner layer location 
where the mean velocity is linear [y + < 5), a location just at the beginning of the 
log-layer region ( y + = 27), and a location in the outer layer at the beginning of the 
wake region (j/ + = 184). The spanwise wavenumber is normalized by the boundary- 
layer thickness 8 C at which the van Driest velocity is As Fig. 2 shows, with 

the exception of the E vv component at y + = 5, spectral amplitude decreases at least 
three decades between the lowest and highest wavenumbers for both the velocity 
and thermal variables. However, the u-energy component exhibits an increase at 
k z 8 c ~ 15 and the overall decrease is closer to two orders of magnitude. The peak in 
the E vv spectrum suggests a periodic r>-velocity component in the spanwise direction 
that characterizes the spanwise vortices imbedded within the boundary layer, fhe 
wavelengths corresponding to the E vv peak at y + = 5 and 27 are AjJ" ~ 100 and 
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50, respectively, 't hese values of Xt are consistent with experimental results (see 
Smits and Dussauge 1996) that show the span wise spacing of streaks in the boundary 
layer to have a mean spacing of A+ = 100. A further validation of the simulation 
can be found in the estimate of the high wavenumber limit of i]/X z (where is the 
Kolmogorov length scale and X z = ‘hrk z ). At y-\- = 27, ?/A zmax is approximately 0.72, 
which can be interpreted as a measure of the relative size of the Kolmogorov scale 
to the spanwise (grid) resolution. While resolution down to the Kolmogorov scales is 
desirable, practical constraints on such a spatial resolution prevents this. 

As Fig. 2(b) shows, only a small spectral region is consistent with the —5/3 
slope. Although other results clearly show a log-law region that equilibrium between 
production and dissipation, the relatively low Reynolds number of this flow precludes 
a large inertial region in the spectrum and precludes a large region where turbulent 
energy and dissipation rate balance. 

Another assessment of the computational accuracy of the simulation is provided 
by the two-point correlation function given by Fq. (4). Figure 3 shows the spanwise 
variation of the correlation coefficients for both the turbulent velocity and thermal 
variables. This distribution can be used to assess whether the computational domain 
in the spanwise direction is sufficiently wide to accommodate the turbulence dynam- 
ics. A measure of the effect is to ensure that the two-point correlation function is 
sufficiently decorrelated over distances L z / 2 (= 0.175 ft 1.63J C ). Figure 3 shows that 
at the y + = 5 and 27 all the turbulent variables have been sufficiently decorrelated 
over distances much less than L z / 2. At y + = 184 the migrations about zero are 
slightly larger, although such migrations are relatively constant for r z > 0.8J C . 

The separation distances at which the minimums in the velocity correlation coeffi- 
cients occur can provide some insight into the compressible turbulent boundary layer 
structure (Kim, Moin and Moser 1987). An estimate of the separation between the 
low- and high-speed streaks can be obtained from the H BU (r 2 )/li u „( 0) autocorrela- 
tion. At y + = 5, shown in Fig. 3(a), R uu (r z )/ R uu ( 0) has a minimum at r z /8 c ft 0.048 
(rf ft 50) which is an estimate of the mean spacing between the high- and low-speed 
fluid at this y-\- location. Similarly, some estimate of the size of any streamwise vor- 
tical structures can be obtained from R vv (r z ) / R vv ( 0). As shown in Fig. 3(a), where 
y + = 5, the minimum in the R vv (r z )/ R vv ( 0) distribution occurs at r z /8 c ft 0.026 
(r+ ft 27). This separation distance then corresponds to an estimate of the mean 
diameter of a streamwise vortex. Interestingly, both these results are consistent with 
the incompressible channel flow results obtained by Kim, Moin and Moser (1987) and 
suggest that, at least in close proximity to the wail, the structure of the compressible 
turbulent boundary is similar to the structure in an incompressible flow. 

3. Results 

File problem under investigation is based on the parameter range of the supersonic 
boundary layer experiment of Shutts, Hartwig and Weiler (1977). In the physical 
experiment, the inlet Mach number was 2.25 and the unit Reynolds number (Re/ in.) 
based on inlet conditions was 6.35 X 10°. Fhe free-stream temperature was 305° R 
and the adiabatic wall temperature was 580° R. 
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Figure 3. Normalized two-point correlation tensor in the homogeneous 2 -direction at 
y + = 5, 27, and 184. 


in the computational experiment, the domain consists of three regions. L he first 
region is fairly well resolved and contains the region of blowing and suction and some 
of the end stage transition phenomena, 2.54 X 10 b < Re x < 4.42 x 10 6 (where Re is 
based on values at the boundary layer edge). 1 he second region also contains some of 
the end stage transition phenomena as well as the turbulent portion of the flow field, 
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4.42 X 10 6 < Re x < 5.67 X 10 s , and has the densest grid in the computational domain. 
The third region is discretized by using a grid that gradually becomes very coarse in 
the streamwise direction, 5.67 X 10 6 < Re x < 9.46 X 10 6 . The total number of grid 
points used in the computation is approximately 47 million (1701 X 55 X 501 grid). 
Unless otherwise stated, the results presented here were determined at the x-station 
corresponding to Re x = 5.548 X 10 6 ( x = 8.8). Based on the wall shear velocity at this 
location, the resolution in the resolved (second) turbulent region is Ax + = 13.9 and 
Az + = 6.8 and the first grid point ofF the wall is at y + = 0.97. The boundary layer 
thickness (based on 99.5 percent of the free-stream velocity) at Re x = 5.548 X 10 6 is 
approximately 1.03 x 10° wall units. 

Figure 4 shows the computed skin-friction distribution compared with a theoretical 
estimate for the fully turbulent regime given by (see White 1974, p. 644 for details) 



where T is the mean temperature, p is the viscosity of the fluid, Re x is the Reynolds 
number based on free-stream conditions, and by 



The region of blowing and suction can be observed in the region 2.8 x 10 6 < Re x < 
3.1 x 10 6 as small variations in the computed skin-friction. Further downstream, 
in the end stage transition region, the skin-friction increases dramatically, reaches 
a maximum and then begins to decrease. As shown in Fig. 4, this streamwise 
development (slope) is coincident with the results predicted by Eq. (6). 

Figure 5 shows the computed van Driest velocity 



normalized by the wall-shear velocity and plotted as a function of the wall normal 
distance (in wall units) at the streamwise location x = 8.8. This location is at a 
distance of approximately 11 boundary layer thicknesses downstream of the end of 
transition. Also plotted in Fig. 5 are the curves U+ = y + and the log-law U+ = 
k - 1 ln(j/ + ) + 5.1 (k = 0.41). The simulation results agree very well with both the 
asymptotic behavior near the wall and the log-law results. 

t he corresponding variation across the boundary layer of the mean density, pres- 
sure, and temperature variables is shown in Fig. 6. Both the mean pressure p and 
total temperature / \ are constant in the inner layer of the boundary layer, but 7 \ 
increases slightly near the outer edge of the boundary layer, the mean density p 
increases across the layer by nearly a factor of two, and the mean temperature T 
decreases by approximately a factor of two across the layer. 
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Re x 10 

X 

Figure 4. Streamwise variation of skin-friction coefficient: Theory from Eq. (6). 



Figure 5. Van Driest mean velocity distribution across boundary layer. Solid line 
DNS, and dashed lines represent U+ = y + and U+ = k 1 ln(?/ + ) + 5.1 in the near- wall 
and log-layer regions, respectively. 







Figure 6. Variation of mean thermodynamic variables across boundary layer. 


With the behavior of the mean velocity and thermal fields across the boundary 
layer separately established, the question arises whether their coupled behavior is 
being properly represented. The coupling of the mean thermal and velocity fields is 
the basis of the Reynolds analogy postulated over 100 years ago. One such measure 
is the Crocco-Busemann relation (White 1971) which can be written for the current 
simulation as 


T 

T e 


7 ' / [2 

-L aw U 

Y7~ 2c p T e ' 


( 9 ) 


where T e is the mean boundary-layer edge temperature, and the isothermal wall 
temperature boundary condition is 1 w = T aw . Note that the recovery factor r is 
often included in the last term of Eq. (9); however, since this equation is derived 

o /q 

under the assumption that the turbulent Prandtl number is Pr t w 1 (r ~ Pr t ), the 
recovery factor was set to unity. Figure 7 shows the DNS data for the temperature 
ratio T/T e as a function of U/U e across the boundary layer for compariosn with the 
theoretical predictions of Eq. (9). The figure shows very good agreement between 
the two quantities across most of the boundary layer. Only in the outer layer, where 
many of the assumptions used in the derivation of Eq. (9) become invalid, do the 
simulation and theoretical results deviate. 

Figure 8(a) shows the behavior of the scaled Reynolds shear stress across the 
boundary layer, and Fig. 8(b) shows the behavior of the corresponding total stress. 
The scaling used in the figure is motivated by Morkovin’s hypothesis (Morkovin 1964) 
which simply states that the “...essential dynamics of these supersonic shear flows will 
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Figure 7. Mean temperature distribution across boundary layer and comparison with 
the Crocco-Busemann relation (dashed line). 


follow the incompressible form...” This observation was quantified by the proposal 
that 

KPwJ 4 \p w J 4 

“would depend little on, Mach number ...” The Reynolds normal stresses are shown in 





Figure 8. Shear stress distribution across boundary layer: (a) scaled Reynolds shear 
stress; (b) total shear stress. Symbol represents incompressible UNS (Spalart 1988), 
and solid line represents current compressible simulation. 
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Fig. 9. It is useful to compare these results with some incompressible data in order 



Figure 9. Reynolds normal stress distribution across boundary layer. Symbols repre- 
sents incompressible DNS (Spalart 1988) at Reg = 1410. Solid and dashed lines are 
current compressible simulation. 

to assess the applicability of Morkovin’s hypothesis and to help validate the present 
DNS data. At the Re x = 5.548 X 10 fa ( x = 8.8) streamwise station shown here, the 
free-stream momentum thickness Reynolds number Reg = (U c< ,0/i' t>o ) is £» 4.25 X 10 3 . 
This corresponds to a wail momentum thickness Reynolds number Reg w = (Uoo9/u w ) 
k, 1.343 X 10 3 . Fortunately, incompressible data exists in this Reynolds number range 
and can be used for comparison. Also shown on Figs. 8 and 9 are the incompressible 
simulation data of Spalart (1988) at Reg = 1410. The results from both the shear 
stress and streamwise normal stress from the current simulation are larger than the 
corresponding incompressible results, particularly in the region of peak values. The 
spanwise normal stress components, ~pio 2 j~p w u\ from both simulations are similar: 
however, the transverse components, pv 2 l~p w u 2 also show that the incompressible 
results are less than the corresponding compressible data. 

From the comparison with the incompressible results, the behavior of the shear 
stress, particularly the peak value, is of concern. Although it can be shown in the case 
of fully developed channel flow that the total stress (turbulent plus viscous) expressed 
in wall units cannot exceed unity, the case of the developing boundary layer is less 
clear. Expressed in wall units, the streamwise mean momentum equation for this zero 
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pressure gradient flow can be written as 




Pv 


dx+ 


d 


dy+ I Oy+ 


P \ — 


Pu 


- =- utT- + H- 


l> \ ^ 




a 


( 11 ) 


where a^ y is the mean viscous stress. Note that the mean density and mean dy- 
namic viscosity ratios appear due to the wail unit scaling that uses the wall values 
for ~p w and fl w {V w ). Since the left hand side of Eq. (11) is small but not negligi- 
ble. no definite statements can be made about the distribution of the total stress 
across the boundary layer. (Recall that for incompressible channel flow the total 
stress varies linearly across the channel.) Nevertheless, as Eig. 8(b) shows, the 
incompressible simulation yields a relatively constant total stress level up to a y + 
value ~ 15 and then decays monotonically with increasing distance to the wall. 1'he 
compressible simulation displays the same trend in the inner layer region as the in- 
compressible data; however, in the log-layer region the total stress exceeds unity, 
increasing slightly and then decreases monotonically as the boundary layer edge is 
approached. While Morkovin’s hypothesis may not apply directly to the behavior 
of the total stress distribution across the boundary layer, the results shown in Eig. 
8(b) suggest that the deviation between the incompressible and compressible results 
may be due to inaccuracies in the compressible simulation, further simulations are 
ongoing (Pirozzoli, Grasso and Gatski 2002)with a different numerical algorithm and 
grid structure that will hopefully shed light on this issue. 

Turbulent fluctuations of the thermodynamic variables p, p , 1\ and (total tem- 
perature) T t are not negligible in such boundary layer flows. Figure 10 shows the 
distribution of these variables across the boundary layer. The pressure fluctuation 




Figure 10. Distribution of root mean square fluctuations of thermodynamic variables 
across boundary layer: (a) non-dimensionalized with free-stream values; (b) normal- 
ized with local mean values. 

levels are relatively constant throughout the inner-layer including most of the log- 
layer region (cf. Fig. 5). As distance from the wall increases, the pressure fluctuation 
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levels decline further (from ~ 2 percent of the mean). Since the mean pressure across 
the boundary layer is relatively constant, the free-stream and local mean scaling have 
very few differences. In the near-wali region, the density fluctuations are also rel- 
atively constant and similar to the corresponding behavior displayed by the mean 
density shown in Fig. 6. Farther from the wall, the density fluctuations reach a new, 
higher plateau (Fig. fO(a)) across the log-layer and a portion of the outer layer. In 
contrast, the normalized density fluctuations (Fig. 10(b)) reach a peak near y + = 20 
and then decrease with increasing distance from the wall. 

As expected, the behavior of the temperature and total temperature root mean 
square (rms) fluctuations are comparable in the near- wall region with both quantities 
increasing linearly with y + . The temperature fluctuations reach a peak value near 
y + = 20 and are about 15 percent of the free-stream mean value. With increasing 
distance from the wall these temperature fluctuations decrease. The same trend is 
exhibited with the normalized temperature fluctuations, although the intensity levels 
are decreased due to local mean temperature exceeding unity across the boundary 
layer (see Fig. 6). Since the mean total temperature is relatively constant and 
greater than unity across the boundary layer, both the free-stream and local mean 
scaling yield qualitatively similar results with only the magnitudes decreasing in the 
normalized case. The total temperature fluctuation levels reach % 5 percent of the 
local mean value. 

Compressible turbulent dilatation effects, such as dilatation dissipation and pres- 
sure-dilatation, do not play a role in wall-bounded flows (Sarkar 1995). Previous 
models (see Gatski 1996) of dilatation dissipation, for example, have assumed a ratio 
of £d/s s ~ O(Mf). Figure 11 shows the dependence of Sd/s s on M t across the bound- 
ary layer. Homogeneous forms of both the dilatation dissipation £ d and solenoidal 
dissipation e s have been used and are given by (Sarkar et al. 1991) 


pe d = -yu 2 k k and pi s = yiOiiOi (12) 

where u>i is the fluctuating vorticity vector. As Fig. 11(b) shows, the turbulent Mach 
number reaches a peak value of « 0.25, whereas the ratio £ d /£ s obtained from the 
simulation data at the corresponding M t is only x 10~ 4 . These results show that 
such models are inadequate for wall-bounded flows and that the dilatation dissipation 
does not contribute significantly to total turbulent dissipation rate. 

4. Conclusions 

A direct numerical simulation of a spatially evolving, supersonic turbulent bound- 
ary layer flow has been described. A fourth-order accurate, upwind- biased, implicit, 
finite-difference method with second-order temporal accuracy has been used in the 
simulation. Power spectra and two-point correlation data have been analyzed to help 
assess the numerical accuracy of the computation. With the exception of the u-power 
spectrum component, these results showed a three order of magnitude decay of the 
power spectral amplitudes across the wavenumber range examined. In addition, the 
decorrelation of the spanwise two-point correlations confirmed that the domain size 
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Figure 11. Influence of dilatation dissipation across boundary layer: (a) ratio of dilata- 
tion to solenoidal dissipation as a function of turbulent Mach number; (b) turbulent 
Mach number variation across boundary layer. 

in the ^-direction was sufficiently large and did not have any adverse impact on the 
current simulation. 

Both the skin friction variation in the streamwise direction and the van Driest 
velocity in the fully turbulent regime agreed quite well with the corresponding an- 
alytical models. The slope of the skin-friction decay in the fully turbulent region 
was consistent with previous theory, and both the near-wall asymptotic behavior and 
log-law were well replicated by the van Driest velocity in the fully turbulent region 
as well. 

The mean turbulent thermal field variables were relatively constant in the near- 
wall region. As distance from the wall increased, the mean density and temperature 
fields showed a corresponding increase and decrease, respectively, across the layer. 
The mean pressure and total temperature remained relatively constant across the 
entire layer. 

Unfortunately, the turbulent shear stress and streamwise normal stress component 
obtained from the simulation data appear to have peak values that are too high when 
compared with incompressible simulation data at similar values of Re$. This result, 
and the relatively low supersonic Mach number, indicate that Morkovin’s hypothesis 
should be applicable. 

Root mean square (rms) distributions of the turbulent thermal fluctuations showed 
that the temperature fluctuation levels can reach as high as 15 percent of the free- 
stream mean value in the inner part of the boundary layer; however, when normalized 
with respect to the local mean, such disturbances are about 9 percent of the mean. In 
contrast, rms density fluctuations are found to have maximum levels that are about 7 
percent of the free-stream mean; whereas, when scaled wfith the local mean value these 
fluctuation levels are about 10% of the mean. Total temperature fluctuations display 
the same qualitative behavior as the temperature fluctuations; how'ever, the maximum 
fluctuation levels, when scaled with the free-stream, are about 10 percent, and the 
local mean scaling yields values in the range of 5 percent. Since the mean pressure 
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distribution was constant across the boundary layer, the relative fluctuation levels 
retain the same relative behavior, whether scaled by the free-tream or local mean 
values. In both cases, the maximum fluctuation levies correspond to approximately 
2.5 percent of the mean. 
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